Fermilab-Pub-04-305-E 
October 27** 2004 



Measurement of the Top Quark Mass In 
All-Jet Events 



D0 Collaboration 



Abstract 

We describe a measurement of the mass of the top quark from the 
purely hadronic decay modes of ti pairs using all-jet data produced 
in pp collisions at ^/s = 1.8 TeV at the Fermilab Tevatron Collider. 
The data, which correspond to an integrated luminosity of 110.2 it 
5.8 pb~^, were collected with the D0 detector from 1992 to 1996. We 
find a top quark mass of 178.5 ± 13.7 (stat) it 7.7 (syst) GeV/c^. 
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The mass of the top quark (m^) is a key parameter of the standard model 
(SM). Knowledge of its value is essential for determining quantum corrections 
to the theory and for limiting the predicted range of Higgs boson masses 012] ■ 
In this letter we report a new measurement of rrit by the D0 collaboration 
in the process pp —>■ ti for the case in which the top and antitop quarks each 
decay to a boson and a h quark, followed by the hadronic decay of both 
W bosons. At lowest order in perturbative QCD, this leads to a final state 
of six quark jets, referred to as the all-jets channel of ti production. 

Measurements of mt have been reported by the D0 and CDF collabora- 
tions based on Run I data (1992-1996) from the Fermilab Tevatron Collider, 
with a pp center-of-mass energy of 1.8 TeV. For the all-jets final state, CDF 
has reported a top quark mass of 186 ± 10 (stat) ± 12 (syst) GeV/c^ [3]. 
CDF, combining leptonic and hadronic W boson decay channels, has mea- 
sured an average top quark mass of 175.9 ± 4.8 (stat) ± 4.9 (syst) GeV/c^ 
The D0 average for mt, based on leptonic decay channels of one or both of 
the W bosons, is 179.0 ± 3.5 (stat) ± 3.8 (syst) GeV/c^ [5^. Combining all 
published measurements yields an average value of 178.0 ± 4.3 GeV/c^ for 
the mass of the top quark [H]. 

The all-jets decay channel has the largest branching fraction of all ti decay 
channels (46%) and is the most kinematically constrained final state, since no 
energetic neutrinos are produced 0. If the jets could be correctly associated 
with their original partons, there would be no ambiguity in the analysis. 
However, the association cannot be made unambiguously. The events of 
interest contain six or more high transverse momentum jets, two of which 
originate from h quarks. The dominant background arises from other QCD 
processes that produce six or more jets. D0 measured the it production cross 
section in the all-jets decay channel to be 7.1 ± 2.8 (stat) ± 1.5 (syst) pb, 
assuming a top quark mass of 172.1 GeV/c^ that was previously determined 
by D0 from other decay channels [?]. This cross section corresponds to 
roughly 360 ti jets events produced at D0 during Run I. 

The D0 detector and our methods of triggering, identifying particles, 
and reconstructing events are described elsewhere [HI E] ■ The measurement 
reported here is based on 110.2 ± 5.8 pb~^ of data from Run I of the Fermilab 
Tevatron Collider, at a pp center-of-mass energy of 1.8 TeV. The events 
were selected with a trigger that required at least five jets of cone radius 
TZ = y^(Ar7)M-~(A0)^ = 0.3, where and r] are the azimuthal angle and 
pseudorapidity, respectively. Each jet was required to have transverse energy 
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(Erp) greater than 10 GeV and |?7| < 2.5. The triggered sample contains 
approximately 2 million events, with an estimated signal-to-background ratio 
of about 1 to 7000. Further selection criteria were applied to this sample, 
including criteria to suppress Main Ring noise, to ensure good jet energy 
resolution, and to remove events consistent with light-quark backgrounds. 
Events containing a reconstructed electron or muon outside a jet cone of 
radius TZ = 0.5 were excluded to avoid overlap with other tt decay channels. 
In addition, we required events to contain at least six jets with Et > 10 
GeV and \ri\ < 2.5. The Et of each jet was scaled to give, on average, the 
correct (true) jet energy. For more information on the D0 jet algorithm, see 
Ref. For information on the jet energy scale correction see Ref. pT] 

These criteria led to the selection of a sample containing 165,373 multi- 
jet events with an estimated signal-to-background ratio of about 1 to 1000. 
Since ti events contain a bb pair, whereas such pairs are relatively rare in back- 
ground events, the signal-to-background ratio can be improved by selecting 
events with at least one jet that may have arisen from the fragmentation of a 
b quark that subsequently decayed to a muon, either directly, or indirectly via 
the chain b —>■ c —>■ fj,. The tagging of 6-jets using embedded muons is effec- 
tive at suppressing background relative to signal because 15-20% of ti events 
are so tagged whereas only 2% of the multijet events in the selected sample 
satisfy the 6-tagging requirements. In this analysis, a 6-jet was defined as any 
TZ = 0.5 cone jet, with Et > 10 GeV and \r]\ < 1.0, that contained a muon 
with pt > 4: GeV/c within its cone. The tag requirement reduced the sam- 
ple to 3,043 6-tagged events, with an estimated signal-to-background ratio of 
approximately 1 to 100. More details on b tagging are given in Ref. p.2j . 

To determine the properties of all-jets events from tt production, exten- 
sive Monte Carlo studies were done. The tt signal events were generated 
using the herwig V5.7 [T^ program and propagated through a detailed 
detector simulation, based on geant V3.15 (13|, and reconstructed with 
the standard D0 reconstruction program. We found that the mean of the 
invariant masses of two triplets of jets formed from the six highest-i^y jets, 
M = |^!!!£i±!!!!2 ~j ^ provided a satisfactory discriminant for distinguishing tt 

signal from background ^3]. We chose the two jet triplets to be those that 
minimized the quantity 



where Miy^ = 77.5 GeV/c^ is the mean value of the reconstructed W boson 




(1) 
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mass in the all-jets tt Monte Carlo events processed through the D0 detec- 
tor simulation and reconstruction programs, and m^^, rrit^ and Myj/^, M^y^ 
are the calculated masses of the reconstructed jets that correspond to candi- 
date top quarks and W bosons, respectively, computed from the jet triplets 
and, within each triplet, the jet doublets. The standard deviations, (Jmt 

^ 31 

GeV/c^ and <7mw ~ 21 GeV/c^, are the average root mean square (RMS) val- 
ues of the mass distributions determined using herwig Monte Carlo events 
generated with top quark masses of 140, 180 and 220 GeV/c^. Minimizing 
provides the correct combination of jets in about 40% of the ti Monte 
Carlo events. 

The top quark mass was measured through the best fit of different ad- 
mixtures of signal and background to the observed mass distribution. The 
fitting technique used is similar to that of Ref. ^^1; which takes account of 
the finite size of every sample in the fit. The posterior probability density 
p{mt, o"ti|Data), computed assuming a flat prior in mass and in the tt cross 
section, is calculated for a set of mass values rrit. For each rrit value, the 
posterior probability density, numerically identical (in this case) to the like- 
lihood L, was maximized by varying au to give the "maximized likelihood" , 
Lmaxijnt) as a function of the hypothesized top quark mass, rrit. The "best 
fitted mass", irifit, was taken to be the location of the minimum of the neg- 
ative log-likelihood curve — hiLmax{f^t)- 

The templates for the top quark signal were generated using a Monte 
Carlo simulation of tt events for a discrete set of masses in the range of 110 to 
310 GeV/c^ in 10 GeV/c^ steps, that is, at 21 mass values. The background 
was modeled using untagged events, that is, multijet data that passed all 
selection criteria except those that define the 6-tag. For each untagged event 
z, a weight Wi is calculated, which reflects the probability of tagging that 
event, such that the sum Wi over all untagged events provides an estimate 
of the background, that is, the number of non-tt 6-tagged events within the 
3,043 event sample. The event weight is the sum of the 6-tag rate per jet, 
which is assumed to depend only on the Et and rj of the jet, and on the muon 
detection efficiency. The tag-rate (t/j) is assumed to factorize as follows: 

tR = T{Et, V, R) = N{R)f{ET, R)g{v, R) , (2) 

where 

f{ET,R) = ao + a,E'T/\ (3) 

g{r],R) = po+Pi\r]\'^ (4) 
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and N(R) is an overall normalization constant. The forms of /{Et, R) and 
g{ri, R) were determined empirically. The tag-rate is divided into 5 bins as a 
function of run number (R) according to the 5 major changes in muon system 
efficiency. For each of the 5 run bins, jets were selected with Et >10 GeV 
and |?7| <1.0. Histograms were made of the Et of the tagged jets and the Et 
of the untagged jets in the data, and their ratio (bin by bin) was fit using 
f{ET,R). Similarly, the distributions of the ratio of tagged and untagged 
{rjl histograms were fit using git], R). The tag-rate was normalized to return 
the number of observed tagged events in the data, tests show very good 
agreement between the tagged data and the predicted background. More 
details on this method are given in Ref. [THj . 

Since the jets in tt events tend to be more energetic, have a more isotropic 
momentum fiow, and have larger transverse energies than those in light- 
quark events, we can enrich the event sample further by event discrimination 
based on a suitable set of kinematic variables. For this analysis, we used the 
following eight variables: Et5 x Etg, \riwi x riw2\, v^, A, S, Nf^J , Hts/Ht, 
and Ht/H, where Eti to Etg and Ei to Eq are the transverse energies and 
energies, respectively, of the six jets, ordered in decreasing Et; t]wi and 
rjw2 are the pseudorapidities of the two hypothesized W bosons; \/I is the 
invariant mass of the A^jeis system; A, the aplanarity, is | of the smallest 
eigenvalue of the normalized laboratory-frame momentum tensor J7] of all 
the jets; iS, the sphericity, is | of the sum of the smallest and next-smallest 
eigenvalues of the same tensor; Nf^^ is the number of jets above a given 
Et threshold, over the range 10 GeV to 55 GeV, weighted by the threshold 
[18J; Ht = EjETj] Ht3 = Ht- Eti - ^T2; and H = EjEj, where the 
sums are over all 71 = 0.5 cone jets with |?7| < 2.5 and Et > 10 GeV. 
Figure^shows comparisons of distributions in each of the kinematic variables 
between background and a tt Monte Carlo signal for nit = 180 GeV/c^. The 
distributions of kinematic variables for events with 6-jets are consistent with 
those without 6-jets. 

The above variables were combined into a single discriminant, calculated 
using a neural network (NN) |I9^ with eight inputs, a single hidden layer 
with three nodes, and a single output Dnn- The network was trained and 
tested with independent samples of herwig Monte Carlo ti signal events, 
at top quark masses of 140, 180 and 220 GeV/c^, and untagged events for 
the background, with the target for background set at and at 1 for the 
signal. Roughly equal numbers of training events were used at each mass. 
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These events (in all, 11,423 for signal and 8,143 for background) were used 
only for training the neural network and producing a single set of network 
parameters. These events were not used in the subsequent analysis. 

Figure El shows D^n for the 3,043 event data sample and for background 
normalized to the same number of events. Also shown for comparison is the 
expected tt Monte Carlo signal for irtt = 180 GeV/c^ multiplied by a factor 
of ten. The final event sample used in the fit to the top mass was defined 
by a cutoff in D^n, which was chosen to minimize the uncertainty on the 
extracted top mass. For a given cutoff on D^n, and a given mass value rrit, 
a distribution was composed by adding the background mass distribution to 
the signal mass distribution, with the signal normalized to the theoretical 
cross section [201 • An ensemble of ~ 100 fake mass distributions was created 
by sampling from the combined distribution. The fitting procedure was ap- 
plied to each fake mass distribution to yield a fitted mass mfu. We thereby 
obtained a distribution of fitted masses, characterized by a mean and an 
RMS, for the given Dnn cutoff and the given value of m^. The procedure 
was repeated for different D^n cutoffs and for top quark mass values of 155, 
165, 175, 185 and 195 GeV/c^. We found that the cutoff Dnn > 0.97 mini- 
mizes the RMS in the fitted mass distributions for the five top quark masses 
considered. 

When applied to the 3,043 events, the requirement of Dj^^ > 0.97 reduced 
the dataset to a final sample of 65 events. Figure 121 shows a comparison be- 
tween the observed mass distribution in the data and the sum of background 
and 175 GeV/c^ top quark signal scaled to the observed number of top events 
(see below). The fitting procedure, with 21 mass values, was applied to the 
observed mass distribution to yield a mass estimate, which was corrected for 
a bias ^Hj of 2.6 GeV/c^ using the relationship 

rufu = 0.712 mt + 53.477 GeV/c^ , (5) 

determined from the Monte Carlo studies. The uncertainties (also bias- 
corrected) were defined to be the 68% interval about the minimum in the 
log-likelihood curve. A systematic uncertainty of 5% arises from the discrep- 
ancy between the jet energy scale in Monte Carlo simulations and that in 
data |llj. The fitting of the tag-rate function introduces a normalization 
uncertainty of 14.9 %. The systematic uncertainty also receives a contribu- 
tion from the bin-by-bin uncertainty due to the limited number of untagged 
events used to model the background. The effect of these systematic uncer- 
tainties on the measured top quark mass was obtained by repeating the fits 
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varying the nominal values by their systematic uncertainties. The effect of a 
small signal contribution to the background sample was checked and found 
not to affect the determination of the top quark mass. 

The insert in Figure El shows the negative log-likelihood as a function of 
the top quark mass for six points near the minimum. After bias correction, 
we find a top quark mass of 178.5 ± 13.7 (stat) ± 7.7 (syst) GeV/c^. As a 
consistency check of our measurement, given the measured top quark mass, 
we can estimate the ti production cross section at that mass. The estimated 
signal in the 65-event sample is 16.6 ± 7 events. This corresponds to a total 
cross section of 11 ± 5 pb, which is consistent with the measured D0 ti pro- 
duction cross section of 5.6 ± 1.4 (stat) ± 1.2 (syst) pb for a top quark mass 
of 172.1 GeV/c2 jH]. 

In summary, we have measured the mass of the top quark, using the 
purely hadronic decay modes in tt events, to be 178.5 ± 13.7 (stat) ± 
7.7 (syst) GeV/c^. This is in good agreement with top quark mass measure- 
ments in other decay channels. 
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Figure 1: Histograms of kinematic variables for ti Monte Carlo signal (solid) 
and background (dashed) normalized to the same area. The Monte Carlo 
signal samples were generated with a top quark mass of 180 GeV/c^. The 
variables are described in the text. 




Figure 2: (a) Djy^ is plotted for data and background. Also shown is the 
Dnn expected for a 180 GeV/c^ top quark Monte Carlo signal scaled up by 
a factor of ten. (b) -Datat in finer bins from 0.90 to 1.05, for data, background 
and 180 GeV/c^ top quark Monte Carlo signal scaled up by a factor of five. 
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Figure 3: Data and the sum of background and Monte Carlo signal plotted 
as a function of the mean mass, M. Insert is — In Lmax as a function of the 
top quark mass. 
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